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Standard heat machines (engine, heat pump, refrigerator) are composed of a System (“working 
fluid") coupled to at least two equilibrium baths at different temperatures and periodically driven 
by an external device (piston or rotor) called sometimes work reservoir. The aim of this paper is to 
go beyond this scheme by considering environments which are stationary but cannot be decomposed 
into few baths at thermal equilibrium. Such situations are important, for example in solar cells, 
Chemical machines in biology, various realizations of laser cooling or nanoscopic machines driven by 
laser radiation. We classify non-equilibrium baths depending on their thermodynamic behavior and 
show that the efficiency of heat machines operating under their influences is limited by a generalized 
Carnot bound. 


I. INTRODUCTION 

Quantum Systems are rarely completely isolated from their environment, whose influence, positive or 
negative, should be considered. The theory of open quantum system was developed IMI to achieve this goal 
and in particular open the way to the study of quantum heat machines, such as engines and refrigerators 
[MCE]. Those models generally assume the interaction between a system and one or two environments, in 
thermal equilibrium, thereby termed heat baths. The efficiency of these machines is limited by the Carnot 
bound, requiring at least two baths at different temperatures in order to extract work. 

Nevertheless, there are many examples in nature where the environment is not in thermal equilibrium, 
such as sunlight, continuous laser radiation, biological cells, etc. Our goal is to establish maximum efficiency 
bounds, as well as to determine the output (work power or cooling power) of quantum heat machines 
operating with non-equilibrium baths. Because the bath is not at thermal equilibrium, the second law does 
not demand the presence of a second bath for work extraction. Our aim is to go beyond a simple situation 
when an environment is a collection of independent heat baths with different temperatures which, in principle 
can be also treated as a single non-equilibrium bath. 

We study a micro or mesoscopic externally driven quantum system, the working fluid, coupled to a large 
environment. On the relevant time-scale, which is longer than the scales of driving and of microscopic 
irreversible processes, the basic parameters of the reservoir are constant, hence the reference state of the 
environment is a stationary state. Such non-equlibrium stationary Systems are well known in macroscopic 
thermodynamics and usually described in terms of local equilibrium with space-dependent temperature, 
density, pressure, fluid velocity etc na. 

However, there are important situations where non-equilibrium character of the environment state is not 
related to spatial non-homogeneity but rather to some internal properties of its state. For example, sunlight 
at the Earth surface is a rather homogeneous environment, the shape of its spectrum roughly corresponds 
to the Planck distribution at the Sun surface temperature T s , but the photon density is much lower than 
the equilibrium one, and the absorption in the atmosphere creates “holes” in the spectrum. 

Another example is a laser radiation m in a continuous wave operation mode, which for the idealized 
single mode situation can be treated as a single quantum oscillator or even a classical monochromatic wave, 
while for multimode case with strong phase diffusion it acts on optically active centers like a non-equilibrium 
bath. Biological machines also provide examples of Systems coupled to non-equilibrium baths either of the 
Chemical nature or consisting of photons or different types of excitons. 

In the next Section we briefly review the theory of quantum heat machines operating with equilibrium 
baths and consider the simplest case, where the working fluid is a two-level system (TLS). In Section III the 
non-equilibrium TLS heat machine is analyzed. The notion of “local temperatures” is introduced. Different 
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examples of non-equilibrium baths are given and they are classified according to their effects on the heat 
machine operation. Finally, in Section IV a general theory for non-equilibrium quantum heat machines 
composed of a general quantum working fluid is given and showed that their maximum efficiency is limited 
by a Carnot-like bound. 


II. STANDARD EQUILIBRIUM HEAT MACHINES 

As a first step and reference point we review the main results of quantum heat machines operating in 
contact with two equilibrium baths. The machine is composed of a periodic modulated working fluid, the 
System, which is permanently coupled to the hot (cold) bath at temperature Tu{c)i and follows a continuous 
cyclic evolution [rarrs) . The total Hamiltonian is 

H tot (t) = H s (t) + + (1) 

i=H,C 

where Hsi is the i-bath free Hamiltonian, Bi its interaction operator and S t is a System operator. The System 
Hamiltonian fulfills periodicity condition Hg(t ) = Hg(t+ 2 ^), H being the modulation angular frequency. 

To illustrate our approach, we choose the simplest realization as an example, nevertheless the same analysis 
may be applied to more complex models. We assume the working fluid is a TLS whose frequency is modulated, 

Hs(t) = (uj 0 + Lo(t))^f-. The modulation may be decomposed into a Fourier series e _ *^o " >dt = ^ ^ q e~' lqQjt 
leading to the Floquet expansion of the Lindblad operator na. 

For a detailed derivation we refer the reader to m- In the interaction picture the evolution equation of 
the working fluid density matrix, p, is given by the quantum Markovian master equation of Lindblad-Gorini- 
Kossakowaski-Sudarshan type pina, which for diagonal matrix elements {p gg ( ee ), ground (excited) state 
population) yields the rate equation 


Pee — {PqG ( uj q)pee + PqG ( UJ q )p gg ) , 

q£Z i=H,C 


where ui q = uio + q£l, P q = ||^ 9 || 2 are the harmonic strength satisfying normalization condition Y2 q ^z Pq = 1- 
The bath coupling spectrum is defined as G l (u >) = e ltul (B l (t)B l (0)}dt and measures the interaction 

strength between the TLS and the bath mode u>. If the bath is at thermal equilibrium, it fulfills the Kubo- 
Martin-Schwinger (KMS) condition (we put h = 1 and ks = 1) [181 ITU] : 


G l (-u}) _ _ U / Ti 

G' M 


( 2 ) 


where T* is the bath temperature. If the TLS interacts only with the i-bath and is not being modulated it 
will equilibrate to this temperature, i.e, p ee /p gg = e~ w °^ Ti . 

We are interested in the thermodynamic behavior at steady state (or limit cycle) where any transient effect 
vanishes. The steady state heat currents from the cold and hot bath, respectively, and the power supplied 
by the source of modulation (work reservoir) are given by the expressions 


qe z ~ r 


P = —Jh - Jc , 


( 3 ) 


where w is the steady state population rate 


w = 



J2 q ezJ2i=H,c PqG l ( w q ) 
J2 q ez J2i=H,c PqG l (u q ) 


and the sign convention is that extracted power is negative. 
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As expected from the Second Law, it has been shown mna that work extraction requires at least two 
equilibrium baths at different temperatures. Nevertheless, this does not apply to non-equilibrium baths. In 
the following sections we address this scenario and find the conditions needed for extracting work from a 
single non-equilibrum bath. 


III. HOW TO INTRODUCE THE FORMALISM OF NON-EQUILIBRIUM BATHS 


We consider the case of a TLS coupled to a single stationary non-equilibrium bath jT2] through the 
interaction Hamiltonian Hi nt = S ® B. Stationarity requires the bath state to be diagonal in the bath free 
Hamiltonian basis. Thermal States are just particular cases of this kind of States. Stationarity is essential 
to ensure that the bath two-times auto-correlations depend only on time difference. 

In order to describe non-thermal but stationary baths, we generalize the KMS condition ([2]) , by introducing 
“local temperature” 


- u /t b ( u ) = G(-uj) 
~ G(lo) 


( 4 ) 


where B is the bath interaction operator. Local temperatures depend on B , the frequency and the state of 
the bath. Only for a thermal equilibrium bath Tb(uj) = T independently of B and w. 

In the next section we show some examples of non-equilibrium baths and calculate their effective local 
temperature. 


A. Harmonic oscillator baths 


Baths like electromagnetic radiation or vibration modes of a material are just a collection of independent 
quantum oscillators with quasi-continuous spectrum of frequencies [2T]. Baths coulcl also be composed of 
fermions, e.g., spin-baths, but even then by suitable transformations (e.g. Holstein-Primakoff Hamiltonian) 
can often be approximated by bosonic baths. The bosonic bath free Hamiltonian is given by 

H B = ^2^ k b^b k , [b k ,b'l]= 5 ki . (5) 

k 

In most applications the bath operator that couples to the System is linear in creation and annihilation 
operators 


B = Ylidkbk + g k bt) ■ 

k 


( 6 ) 


The state of the bath is assumed to be stationary and hence diagonal in the particle number basis. Then, 
the coupling spectrum yields 


GM 


\ 9 k\\n k + -lo), w > 0 1 

Y^ k \gk\ 2 n k 5(L0 k - w), io < 0 J ’ 


( 7 ) 


where n k = Tr (pBb k b k ) is the k-mode population, the upper (lower) line in |7| is the emission (absorption) 
rate. The local temperature is given by 


or 


T b (oj) 


ln [ 


u: 


gM i 

G(-w)J 


( 8 ) 


e -a,/r B (o,) = n ^) 

n(w) + 1’ 


(9) 


where n(u>) = 


V,. \gk\ 2 nkS(u k -uj) 


denotes the average population number for the frequency w. 


Y2k 1 9k\ 2 S(u k -u) 

As we show below, the frequency dependence of the local temperature determines whether work can be 
extracted from the single non-equilibrium bath. 
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Figure 1: (Color Online) Local temperature of a TLS coupled to sunlight. The continuous line corresponds to the 
case of thermal radiation at the temperature of Sun surface T s ~ 6000/\, while the dashed line to sunlight out of 
equilibrium due to a geometrical factor, A = 2.5 x ICC 5 (see main text). 


1. Sunlight 

The Sun is a thermal source, emitting thermal radiation at T s = 6000&. Due to geometrical considerations 
122 [23], just a small fraction of the emitted photons reach the Earth, reducing the effective mode population 
and thereby n(w) = A[e“/ Ts — 1] , where A = 2.5 x IO -5 is a geometric factor equal to the angle subtended 

by the Sun seen from the Earth. Effectively, sunlight on Earth is out of equilibrium, and Systems with 
different frequencies, will “equilibrate” to different temperatures. Fig|T[ shows the equilibration temperature, 
Pee/Pgg = as a function of the TLS frequency. Moreover, the atmosphere acts like a filter and 

produces a more complicated shape of n(cu) with many “holes". 


2. Multimode laser radiation 

A multimode laser radiation in a continuous wave operation mode may be modeled as a bath at the phase 
average, ps, of the multimode coherent state ps, 


k 

Pb = t4|vac)(vac|Z7.J, |vac)(vac| = (^) |O fc )(O fc |, (10) 

where U' z is the displacement operator 


Ulb k U z = b k - z(uj k ). (11) 

This non-equilbrium bath is obtained by displacing a thermal equilibrium bath at zero temperature and 
performing a phase averaging. The phase average of the bath state, or diagonality in the photon occupation 
number basis, is required in order for the bath to be stationary. The local temperature for this bath is equal 
to 


co 

ln(l+ |z(w)|- 2 )' 


T b M 


( 12 ) 
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For large |z(w)| yields 

T b (u;)~u;|zM| 2 . (13) 

Therefore, for a constant shift, z, for ali the modes, the local temperature is a linearly increasing function 
of the frequency. 


3. Squeezed thermal bath 

The state of a stationary squeezed thermal bath is the phase average 122125], Pb, of the following density 
matrix 


p B = Z- x S r e- HB ' T ^Sl , 

where S r is the squeezing unitary operator defined by 

Sj.bkS r = cosh (r(uj))bk + sinh(r(w))6^. 
Its local temperature is given by the expression 

T b (u) = 


U! 


(14) 

(15) 

(16) 


ln{l + [n Te9 (w) + (2 n Tcq {u) + 1) sinh 2 (r(w))] } 
where riT eq (co) = [e“/ Te<! — l]” 1 , is the mode population without squeezing. For large r* (161 reduces to 

T B ( W )~|(2n r .>) + l)e 2 V (17) 


B. Classification of non equilibrium baths 

The functional dependence of T B {u>) can be used to classify non-equilibrium baths. It depends on the 
bath state, coupling to the system and frequency. In this sections we consider a heat engine composed of 
a TLS coupled to a single non-equilibrium bath and find the required conditions for work extraction. The 
dynamics of the total system is governed by the Hamiltonian 

Htot.it) = (w 0 + ut it)) ——b S <S> B + H b ■ (18) 

The steady heat current is a sum of contributions corresponding to all harmonics of f l (see section IV for a 
general case) 


J =E 


UqPq 

W + 1 

qkz/L 

with the TLS steady state population ratio 


(G(-UJ q ) - G(u q )w ) = -P, 


w = 


Pee 

Pgg 


ss 


Yhq&Pq^i Uq) 
Sqe z PqG{u q ) 


(19) 


( 20 ) 


1. Non-equilibrium but passive 

Using Eqs[4 ( 191, (p0| we can write the expression for power supplied to the system as 


P = z 


1 E (?i-®)nP, 1 P, 2 GK)G( W , 2 )(e-“-^K) 


_ p ^qi /Tb (bJqi) 


{qi>q2t%} 


)• 


( 21 ) 








where 
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^ = ^2 P q G(uj q ) [1 + e-^/^K)]. 

q£ Z 


( 22 ) 


A sufficient condition 


w 92 /Tb(w 92 ) > U} qi /T B (u} qi ), for {q 1 >q 2 } 1 


(23) 


assures that P < 0 and hence the engine extracts work from the bath. 

We define the passivity function as /(w) = ( u/Tb(u !)). If for all frequencies /(w) > 0, no work can 
be extractecl. We term such couplings to bath passive (in analogy to passive States [2DJ E5H2B] which do 
not allow for work extraction). Previous examples, based on linear coupling to bosonic bath, are passive if 
we consider only frequency independent deformations of a thermal bath (constant filtering, displacement or 
squeezing). From Eq 21 we deduce that work extracted from a single non-equilibrium bath depends on the 


coupling spectrum shape and requires the passivity function /(w) to be negative in some range of frequencies. 


2. Two equilibrium baths as a single non-equilibrium bath 


As a first example of a non passive bath, we consider the Standard quantum heat engine where the working 
fluid interacts with two baths at equilibrium. We model it as a quantum heat machine with a single non- 
equilibrium bath. For this we need to consider the spectrum as a sum of two baths spectra. Therefore, the 
local temperature of the composed bath satisfies the following relation 


e -u>/T B (u) 


e -“/ Th G h (pj) 
G h (uj) + G c (uj) 


C-u/TcQC^ 

G h (u;)+G c {u>) 


= e~“ /Th {l - m(w)) + e - “/ Tc m(w), 


0 < m(w) < 1, 


where we use the fact that both baths are in equilibrium and the Standard KMS condition holds. The effective 
Boltzmann factor is a weighted average of both two Boltzmann ones, where the weights depend on how strong 
is the working fluid coupled to each bath at the given frequency. Therefore, T > T n (ui) > T c . For the 
sake of simplicity, we assume that the bath coupling spectrum overlaps with only two harmonic frequencies 
(u> qi > uj q2 ). As shown in m this condition is required in order to achieve high efficiency. Work extraction 
requires the hot bath being coupled more strongly to the high frequency mode ( G h (uj qi ) G c {uj qi )) and 
the opposite for the low frequency mode ( G h (ui q2 ) <C G c (oj q2 )) P3|. The efficiency of the engine is given by 
the Carnot-type formula 


T B{u q2 ) 
Tb{u q J 



The extreme case, where the engine reaches the maximum efficiency, the Carnot bound, is when the bath 
are spectrally separated, Tb( ui qi ) — Th and Tb(cu 92 ) ~ T c . 


3. Non-equilibrium and non-passive bosonic bath 


As shown above, frequency independent deformation of a thermal bosonic bath with linear coupling, creates 
passive baths. Therefore, work extraction requires the use of selective filters. In order to not contradict the 
second law of thermodynamics, which forbids work extraction from a single thermal bath, this selective filter 
should involve the presence of other bath, a non equilibrium process or a hidden work injection. 

Assume a single bosonic bath at the equilibrium temperature T eq and linearly coupled to the System. If a 
selective filter, A(w), is applied, the local temperature satisfies 


■u/t b (u) = 

A (ui)n(ui) + 1 


(24) 


What are the conditions required for this filter to allow work extraction? 
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Consider again the bath coupling spectrum which overlaps with only two harmonic frequencies u > qi , oj, 


with qi >^ 2 - As shown in m this condition is required in order to achieve high efficiency and due to ( |21[ ) 
the power is given by 


J Qr> 


P = z 1 {qi~ q2)ttP qi P q2 G(uj l 


Q i 


)G(uj q2 ) 


,JTb (a 


2 ) _ g ^<21 fT]3 (^qi ) 


)• 


(25) 


Assume that we reduce the population of the mode oj Q2 (A(ol> 92 )<1), and we do not filter the other mode 
(A(w gi ) = !)• In order to allow work extraction, the filter should satisfy the following condition 


< 1. 


As paradoxical it may sound, by reducing a specific mode population, we can extract work from a single 
thermal bath! The filtering lowers the effective temperature Tb(w 9 ,), reducing the excitations that are 
emitted to this mode, an “saving” energy which is ultimately transformed into work. The efficiency of this 
machine is bounded by a generalized Carnot limit 


V = 



T B{Uq 2 ) 

T eq 


where T eq and Tb(w 92 ) play the role of the effective hot and cold bath temperature respectively. 


4. Deviation from equilibrium of engineered bath 

Using the selective filtering introduced in the last section, we can engineer a non-equilibrium bath from 
an equilibrium one, and characterize its deviation from equilibrium by the parameter 

D= 1 - A (ujq 2 ). 

For the equilibrium bath D = 0, and when we start reducing population of the modes with frequencies 
around w 92 , the bath will go away from equilibrium producing some non-equilibrium effects, like frequency 
dependent equilibration temperature. Nevertheless work extraction will be possible only when the bath is 
far enough from equilibrium, i.e. 

D > n(ui qi )(e ulqi ^ Tqq - e ulq2 / T ' !q ). (26) 

When the deviation from equilibrium increases, the local temperature of the lower frequency mode reduces 
and the efficiency of the quantum heat engine rises. For D < 0 the bath is also out of equilibrium. But, in 
this case instead of reducing the mode population, it is being increased by some external mechanism (for 
example, selective concentration of light). The equilibration temperature of the system will depend on the 
frequency, but for D < 0, work cannot be extracted. Again, paradoxically, the increase of energy in “incorrect 
modes” reduces the possibility of work extraction. The bath is taken away from thermal equilibrium, but in 
the “opposite direction” to that leading to work extraction. 


IV. GENERAL THEORY FOR NON-EQUILIBRIUM QUANTUM HEAT MACHINES 

The model based on the TLS and studied above is an example of a large class of open quantum Systems with 
physical Hamiltonian Hg(t) = Hg(t + t) under the assumption that the perturbation frequency fl = 27 t/t 
is comparable to or higher than the relevant Bohr frequencies. Such fast driving is typically providecl by 
a strong coherent laser field and appears in the thermodynamical approach to the theory of lasers [131223- 
EU, various types of laser cooling [32] Ea, optomechanical clevices EU, etc. Potential applications include 
new light harvesting Systems, both of biological nature or man-made devices. For this class of Systems a 
consistent theory can be developed which includes the general case of stationary nonequilibrium environment 
characterized by local temperatures. The laws of thermodynamics can be derived and Carnot-type bounds 
are obtained. 






A. Master equations 


We begin with a presentation of the canonical construction of the Markovian generator for an open system 
weakly coupled to a stationary, but generally nonequilibrium, environment. The system is assumed to be 
“small” and described by the periodic in time plrysical Hamiltonian Hg(t) = Hg(t + r) under the assumption 
that the perturbation frequency f 1 = r/r is comparable to or higher than the relevant temporal Bohr 
frequencies. We assume that the Hamiltonian of the system already contains all Lamb-like shifts induced by 
interaction with environment na amni- 
The system bath-interaction is parametrized as 

Hint = ]T ® .Bo, (27) 


where S a and B a are hermitian operators of system and bath, respectively. The environment (bath) is 
assumed to be a large quantum system with practically continuous Hamiltonian spectrum and a proper 
behavior of multi-time correlation functions of relevant observables. The initial state of the bath is stationary, 
therefore is invariant with respect to the free dynamics of the bath and satisfies ( B a )g = 0, where ■ - )b 
denotes the average over the bath state. We assume also, for simplicity, that the cross-correlations between 
B a and Bp vanish for a/,0. 

Applying the Floquet theory one obtains the following decomposition of the associated unitary propagator 


U{t) = Texp |—i J 77s(s)dsj = P(t)e 


—iHt 


(28) 


where P(t) 


P(t + r) is a family of periodic unitaries and H is the averaged Hamiltonian satisfying 


U (t) = e~ iHT . 

The Floquet operator U(t), and the averaged Hamiltonian H posses common eigenvectors {</>&}, i.e., 

H(j) k = e fc </> fc , U{r)(j) k = e~ ltkT <j> k , 


(29) 

(30) 


where {e^} are called quasi-energies of the system. These properties imply a particular form of the Fourier 
decomposition 


where 


S a (t) = U(t)'S a U{t ) = ^ S a (Q q )e- m ', 


(31) 


{ w ?} — {w + qtt} ; {w} — {e k — e;}, q £ Z}, (32) 

i.e., it is a set of all sums of the relevant Bohr quasi-frequencies and all multiplicities of the modulation 
frequency. By {u> g }+ we denote the subset of {w 9 } with non-negative relevant Bohr quasi-frequencies u. 
The operators S a (ui q ) are subject to relations 

B a (ca<j)^ — S a ( ^q); 

[■ H,S a (uj q )] = -wS a (u> q ). (33) 

Physically, the harmonics qfl correspond to the energy quanta which are exchanged with the external periodic 
driving. 

Repeating the construction of the weak coupling generator in the interaction picture one obtains 

( 34 ) 

“ {“ 5 }+ 

with a single term defined as 

=^{Ga(w 9 )([»S' Q (w g ),pS ct (a; 9 ) t ] + [S a (u> q )p, S a (u! q )^ 

+ G a (-u q ) ([S a (u q )\ pS a {u q )\ + [S a (u] q )^p, S a (Q q )]j |. 


(35) 
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Using (33 ( one can show that C ^ commutes with —i[H, ■] and possesses a Gibbs-like stationary state 


q /ui)H/T a (ui q ) 


(36) 


with the local temperature T a {uj q ) corresponding to the coupling channel (a,u g ). The “renormalizing" 
term u> q /u> in front of H in the Gibbs-like state takes into account the total energy exchange including 
qfl corresponding to external clriving device. The properties (33 ( imply that the generator (35l transforms 
independently the diagonal and off-diagonal elements of p computed in the eigenbasis of H. 

The MME the Schroedinger picture possesses the following structure 


dpsch (t) 

dt 


= -i[H s (t),p ac h(t)} + (U(t)C U(ty)p sch (t ) 


where U(t)p = U(t)pU(ty. A very useful factorization property for the solution of (37) holds 

Psch{t) = U(t)e tc p sch {fi), 


(37) 


(38) 


which allows to discuss separately the decoherence/dissipation effects described by C and the unitary evolu- 
tion U{t). 


B. The Laws of Thermodynamics 

The structure of MME’s derived above and the introduced notion of local temperatures allow to formulate 
the first and second law of thermodynamics in terms of energy, work, heat and entropy balance. The basic 
tool is the following inequality valid for any LGKS generator C with a stationary state p [5] and arbitrary p 

Tr(£p[ln p — lnp)]) < 0. (39) 

This inequality is used to show positivity of entropy production under the assumption that the physical 
entropy of the System S(t) is identified with the von Neumann entropy of its state 

S(t) = -Tr(p(t)lnp(t)). (40) 


1. Entropy balance and local heat currents 


In the case of fast driving there is no obvious definition of the temporal internal energy of the System 
because a fast exchange of energy quanta qVL between system and the external source of driving makes a 
temporal partition of energy between both Systems ambiguous. The situation is different for the entropy 
balance because the entropy change is due to irreversible processes which are slow under the weak system- 
environment coupling assumption. The weak coupling scheme yields the coarse-grained in time effective 
dynamics described by the MME (3TI what suggest the following definition of the heat current J Q ( uj q ) 
supplied to the system by the coupling channel (a,uj q ) and involving the averaged Hamiltonian (multiplied 
by uj q /dj) 


J° q (t) = ^Tr (HC? qP (t)), 


(41) 


where p(t) is the system density matrix in the interaction picture and according to (371 given by 

p{t) = e tc p{ 0). (42) 


Those definition allow to formulate the Second Law which is again a consequence of (391 applied to each 
single coupling channel 


dt 


{w,}_ 




(43) 
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where <Jq (t) > 0 is an entropy production caused by a single coupling channel (a, oj q ) given by 

<(i) = Tr(C% q p{t)[\np{t) - lnpgj) > 0. (44) 


For the case of environment composed of several independent heat baths the equation (43) reduces to the 
Standard form of the Second Law for open Systems with usual temperatures. 


2. Steady state regime 


Under natural ergodic conditions and due to (371 any initial state tends to a limit cycle (or fixed point in 
particular cases), i.e., 


p(t) -» p(t) = U(t)pU(t)^ = p(t + t), where Cp = 0. 
Then the entropy S{p(t)) and the heat currents given by 

= |Tr (HC? q p), 

become constants leading to the following form of the Second Law 

T9 

,,, 


E E 


{“<jU 


T a (iOq) 


< 0. 


(45) 


(46) 


(47) 


The averaged internal energy of the System is constant in the limit cycle, and hence we can use the total 
energy conservation to write the First Law in the form 


P = -Y. E 

“ {“?}+ 


(48) 


where P is the stationary power and if it is negative, it is supplied to the source of external driving. 

Remark 

Because each Cq transforms diagonal (in H basis) elements of the clensity matrix into diagonal ones, the 
stationary state p is diagonal and hence the expressions for the stationary local heat currents and power 
involve only diagonal elements and the “classica!” transition probabilities between them. 


C. Carnot bound at steady state 


In the steady state regime the incoming and outgoing heat currents can be defined as follows 

j (+) = E J :.< J<_, = E [-<!■ < 49 ) 

{a,{w ? } + ;J“ a >0} {a,{w 5 } + ;J“ 9 <0} 

We can introduce also effective “hot/cold bath temperatures” by averaging the inverse local temperatures 
with the weights proportional to incoming/outgoing heat currents 


(50) 


1 V 

J OJ q 

1 1 V 

JS 

u q 

1 

T(+) ^ 

{“,{w,}+;J“ 9 > 0} 

UMi 

r Q (w,)’ r(-) 

{a,{u q }+;J2 q < 0} 

L Jc-) J 

T a {u) q ) 


Combining now (471 - (501 and the Standard notion of an efficiency of a heat engine one obtains the generalized 
Carnot bound 


_p r(-) 

V ~ J(+) - _ T(+) ’ 


(51) 


which again coincides with the Standard one in the case of environment composed of two heat baths. 
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V. CONCLUSIONS 

We showed that quantum machines weakly coupled to a single non-equilibrium stationary environment, 
and subject to fast periodic driving by work reservoirs, can be described by the thermodynamical principies 
and bounds which are very similar to the Standard ones if only the proper definitions of the basic notions 
are used. In particular the notion of local temperature which depends not only on the state of environment 
but also on the form of system-environment coupling is crucial. 

The developed non-equilibrium theoretical framework may be used also to described the Standard heat 
engine model, which operates under the interaction with two thermal baths. They can been effectively 
described as a single non-equilibrium bath. Therefore, we show that Standard heat engines are just particular 
examples of non-equilibrium heat machines. 

Starting from a bosonic thermal bath we showed how to obtain non-equilibrium baths by using different 
filters. Displacement or squezeeing operations may also be used. We found out that such non-equilibrium 
baths may be divided into two different types: (i )passive, which equilibrate Systems to different temperatures 
depending on their frequency but cannot drive heat engines, such baths can be obtained by frequency inde- 
pendent transformations of equilibrium ones; (ii) non-passive , which in addition allow work extraction from a 
single bath and, in the case of bosonic reservoirs, can be engineered by frequency dependent transformations 
of equilibrium States. They are also farther away from equilibrium than passive baths. 

A case where non-equilibrium bath is not stationary but, for example, is also perturbed by an external 
periodic driving is another interesting topic with possible applications. A natural example is a spin-1/2 
coupled to a spin-bath, both periodically perturbed by external magnetic field. It seems that the theory 
presented above can be extended to these cases as well. 
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